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1. Introduction 



Sh(s) 



Understanding the complex behaviour of glass-forming systems cooled to below their 
glass transition temperature remains a challenge. At first sight, the problem might 
appear to be hopelessly difficult, since time-translation invariance is in general broken 
and the properties of observables may hence depend on the prehistory of the material 
under study (i.e. thermal, mechanical,. . . ). On the other hand, an important discovery 
has been the observation of dynamical scaling, see [1], which occurs quite independently 
of whether the equilibrium state is critical or not. In recent years, it has been realized 
that many aspects of dynamical scaling are conveniently first studied in non-disordered, 
i.e. ferromagnetic systems. After a quench to or below their critical temperature T c , 
these systems undergo an ageing behaviour which in many respects is quite similar to the 
one in glassy or kinetically constrained systems. The manifold problems which arise in 
the study of ageing in simple magnets or glasses are reviewed, e.g. in [2, 3, 4, 5, 6, 7, 8, 9]. 

For notational simplicity, we shall in what follows consider magnetic spin systems. 
Convenient tools for the study of ageing behaviour are the two-time autocorrelation and 
autoresponse functions 

C(t,s) = (<f>(t)cf>(s)) ~s- b f c (t/s) (1) 

~ s-^fRit/s) (2) 

h=0 

where 0(t) is the order parameter at time t and h(s) is the conjugate magnetic field 
at time s. The scaling behaviour is expected to apply in the so-called ageing regime 
where t, s ^> t m i cro and t — s ^> t m i cro , where t micro is a microscopic time scale. 
Furthermore, we tacitly assumed that the scaling derives from the time-dependence 
of a single characteristic length-scale L(t) which measures the linear size of correlated 
clusters. In this paper, we shall always consider algebraic growth laws, viz. L{t) ~ t 1 ^, 
where z is the dynamic exponent. Then the above forms define the nonequilibrium 
exponents a and b and the scaling functions fc{y) and fn(y)- For large arguments 
y — > oo, one generically expects 

fc(y) ~ y~ Xc/z , f R (y) ~ y~ XR,z (3) 

where Ac and A#, respectively, are known as autocorrelation [10, 11] and autoresponse 
exponents [12]. This description applies to many simple magnets quenched to a 
temperature T < T c as is well-known, see [2, 3, 5].t On the other hand, for glasses 
quenched to below the glass-transition temperature, a slow cross-over between an 
algebraic growth law at short times towards a slower (logarithmic) growth at larger times 
appears to give a better description of the data, see [8] and references therein. Recently, 
evidence was found that for spin glasses quenched onto their critical temperature T = T c 
a simple scaling of the two-time observables compatible with an algebraic growth law 
applies [13]. 

f Exceptions occur for example in the 2D XY model with a fully disordered initial state. 



Scaling and ageing of the alternating susceptibility in spin glasses 3 

In recent years, from the study of ageing in simple magnets it has been proposed 
that global dynamical scaling might be extended to a local scale-invariance [14, 15]. 
One of the simplest predictions of that theory is the explicit form of the two-time 
autoresponse function. It readsf 

R(t, s) = s-^fnit, s) , f R {y) = / y 1+af -^'(y - (4) 

where a' is a new independent exponent and /o is a normalization constant. The 
independent existence of the exponent a' was recognized recently [17] for z = 2 and 
we shall extend that argument to arbitrary z in appendix A. Previous derivations of 
fn(y) [14, 15] had assumed a = a' from the outset. The only known example with 
distinct exponents a ^ a' was the ID Glauber-Ising model at temperature T = and 
with initial power-law correlations of the form (crjUj) ~ \i — j\~ v with v > 0. The 
exact solution for R(t, s) [18, 19, 20] is of the form (4) and one can read off a = 0, 
a' = — \ and Xr = 1 [17]. In appendix B, we give a second example and show that the 
treatment of phase-ordering in d spatial dimensions at zero temperature in the spirit of 
the OJK-approximation [21, 22] leads to 

_d-\ , _d-2 _d 
aojK - 2 , a OJK - - , A^ojk - ^ 

The same scaling function and hence (5) is also found from the gaussian approximation 
to phase-ordering kinetics [23]. § 

Eq. (4) with a = a' is perfectly reproduced in simple magnets undergoing coarsening 
after a quench to T < T c , most notably the 2D 1 3D Glauber-Ising model [15, 24], the 
3D XY model [25] and in several exactly solvable systems [14, 27, 17]. || For critical 
quenches T = T c , the agreement between numerical data and (4) with a = a' is almost 
perfect in the 2D/3D Ising and XY models [15, 27, 25, 28], the ID contact process 
[29] and is exact in several exactly solvable models with z = 2 [14, 12, 17]. On the 
other hand, a second-order e-expansion calculation from renormalized field-theory gives 
a small but systematic deviation with respect to eq. (4) [9, 30]. 

Of course it would be interesting to see whether a scaling description or even an 
extension to local scale-invariance might be applicable to glassy systems as well. Such 
a comparison may be far from straightforward, however. Indeed, a possible scaling 
behaviour of two-time correlators C(t, s) and integrated responses J du R(t, u) have been 
discussed since a long time, both theoretically and experimentally, see [1, 3, 6, 4, 31, 8] 
and debates have arisen over several central issues of which we mention a few. First, 

| Technically, this requires that the order parameter is a quasi-primary field under local scale- 
transformations [14]. This concept is the analogue of the one used in conformal field-theory [16]. 
§ We observe that in both examples a = a' + i but it is still open to what extent this might be a 
general relationship. 

|| A further extension of local scale-invariance with z = 2 in d spatial dimensions to a new type of 
conformal invariance in d + 2 dimensions yields a prediction of C(t,s) which is in agreement with 
numerical data in the 2D Glauber-Ising model [26] . 



(5) 
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for glasses quenched to below their glass temperature, it is not entirely clear whether 
the growth law for L(t) is algebraic or logarithmic [32]. Second, even if an asymptotic 
power-law scaling is accepted, there has been an intense debate on whether scaling 
occurs according to the so-called 'full ageing' scenario, that is in terms of the scaling 
variable y = t/s (possibly with small logarithmic corrections) [33], or else if a 'sub- 
ageing' scenario applies, with a scaling variable £ := [t 1 "^ — s 1 ~ M ]/(l — /x), where /i is a 
free parameter [31]. The usual power-law scaling is recovered in the /i — > 1 limit, but in 
many experiments the data are fitted with values of /j, as low as 0.8 — 0.9. It has been 
suggested recently [34] that values of /J, < 1 merely result from a quench to below the 
glass transition which is not yet sufficiently rapid, but the repetition of that experiment 
on other substances has not yet led to unambiguous conclusions [31]. Third, it is even no 
longer obvious that the commonly studied spin glass models really mimic sufficiently well 
the experimentally studied materials (in spite of well-established qualitative similarities 
[35]): recent simulations on 3D/4D Ising and Heisenberg spin glasses provide evidence 
for cumulative ageing and rejuvenation phenomena in temperature cycling which are 
not observed in real spin glass materials [36]. 

In view of these many difficulties, it might be simpler to consider the behaviour of 
glassy systems from a different point of view. One rather works with a time-dependent 
(oscillating) magnetic field and studies simultaneously the dependence on time and 
on the imposed oscillation angular frequency u>. For a harmonic magnetic field, it is 
common to consider the real and the imaginary part of the magnetic susceptibility 



where R(t,s) is the linear response discussed above. In this setting, 1/uj provides the 
second time-scale, the natural scaling variable is y = out and the scaling regime should 
be reached in the limit u — > and t — > oo. In many experiments and simulations, one 
averages over at least one period of the oscillating field, see e.g. [35]. Then in a great 
variety of glass-forming substances quenched to below or near to their glass transition 
point one observes good but not always perfect evidence for an wt-scaling behaviour of 
the following form of the period-averaged dissipative (imaginary) part [32, 37] 



where x'L * s thought of as a 'stationary' contribution while the ageing behaviour is 
described by Xage- The amplitude A' a ' gc and the exponent b" are obtained from fits to the 
experimental data but there does not seem to exist a relationship with the exponents 
a, a', b, \c,r in the literature. 

Similar scaling forms have been proposed for the dispersive (real) part x' but in 
practice the imaginary part is usually easier to measure. It is usually thought that 




(6) 





b' 



b". 
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In this paper, we shall try and see whether a relationship between the exponents 
b" and b' of the alternating susceptibilities (6) and the other nonequilibrium exponents 
arising in the two-time observables can be found. Assuming the validity of local scale- 
invariance for the relatively large values of z (usually, z ~ 5 — 7 [8]) found in many 
studies of spin glasses and hence the explicit form (4) for the autoresponse function,^ 
we shall show that for large times 

x >, t) = x » + ^Y 2 'M + o {t- XR,z ) 

X V, t) = X '» + + o (t' x ^ z ) (8) 

and we shall calculate the scaling functions x'i an d X2 explicitly in section 2. For the 
asymptotics of these scaling functions we expect X2(y) ~ V~ b ' an d x'iiv) ~ V~ b " f° r 
y ^> 1 and obtain the relation 

b' = b" = a - a' (9) 

This should be compared to the experimentally found scaling (7) of XageO^) anc ^ 
similarly for x' agc (cut) . In section 3 we compare these results with Monte-Carlo data 
from the 3D Ising spin glass with a binary distribution of the couplings and discuss to 
what extent the scaling relations (9) and the explicit scaling functions agree with existing 
experimental data. We conclude in section 4. In appendix A we derive eq. (4) from (an 
extension of) local scale-invariance and in appendix B we revisit the OJK-approximation 
of coarsening kinetics and derive (5). 

2. Scaling of the alternating susceptibility 

We now analyze the frequency-dependent scaling of the alternating susceptibility. For 
notational simplicity, we concentrate first on x"(w,t) as given by eq. (6). In order to 
make the scaling behaviour explicit, we must convert this into a more convenient form 
which can be done as follows [38]. We observe that the time difference r = t — u plays 
a central role since depending on its value either an equilibrium behaviour or else an 
ageing behaviour is obtained. Specifically, it can be shown [39] that there is a time-scale 
t p ~ ^ with < ( < 1 on which the transition between the two regimes occurs such 
that R(t, s) ~ R cq (t — s) for t — s < t p and R(t, s) = s~ 1 ~ a f R (t/s) as given in eq. (4) for 
t — s > t p . + In addition, one measures for u ~ t the response with respect to a change 
in the initial conditions and then instead of (4) one expects R xs R- m i(t) ~ t~ XR l z [2]. 
We therefore must introduce a further time-scale t £ such that t — 1 £ = 0(1). Changing 
variables and then splitting the integral into three terms corresponding to these three 

Our recent study of the critical 3D/4D Ising spin glass [13] showed that the form of the measured 
thermoremanent magnetization p(t, s) = J^du R(t, u) agrees with the prediction of local scale-invariance 
for t/s < 20, which is also in the sector encountered in wt-scaling. 
+ Explicitly, ( = A/(d + 2) in the d-dimcnsional spherical model [39]. 
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regimes, we have 

x"(uJ,t) = / dr R(t, t — t) sin lot 
Jo 

= I dr R(t, t — r) sincjr + / dr R(t, t — r) sin cur 
Jo Jt p 

+ / dr R(t, t — r) sin lot 

Jt e 

+ f Afl/ * / drc sincjr 

= X -;V) + r- j\ h (jl-) + o (10) 

In the third line, we used the asymptotic forms of R(t, s) as described above. This 
means that the cross-over between the equilibrium and the ageing regimes is assumed 
to be very rapid. In the last line, we restricted ourselves to the long-time limit t — > oo. 
Here, the function Xii 00 ) merely depends on the equilibrium form of the linear response 

Req(t,s). 

In this way (analogously for x!) the scaling form (8) is obtained. This derivation 
also shows that the often-found stationary term in the integrated response [3, 40, 4, 31, 
32, 8, 41, 34, 37] does not require the separation of a similar 'stationary' part in the 
response function R(t, s) itself. 

We now analyze the second term in the above expression for x" ■ Using the explicit 
form (4) for the scaling function f R , we have 

x >, t) = X » + U' a S + O (t~^') (11) 

where expansion of the sine followed by term-wise integration gives 

S := [ dv (l-v)~ 1 ~ a+XRlz v~ 1 - a ' sm(uutv) 
Jo 

^(2n+l)! V ' * J ^ ] 

= B (l-a',— -a) ut (12) 



z 



1 — a' 2 — a' 3 1 — a — a' \ R 2 — a — a' X R oj 2 t 2 
x 2-F3 I — - — , — - — ; - \- 



2 ' 2 '2' 2 2z' 2 2z' 4 

where 2F3 is a hypergeometric function, 

B(z,w) = ffiH^ = C&uu'-^-uY 3 - 1 (13) 



+ w) Jo 
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is Euler's betafunction and the identity T(2z) = 2 2z 1 T(z)T(z + 1/2)/ y/n was also used. 

We proceed to extract the leading behaviour for large values of the scaling variable 
y = ut. Recall first the asymptotic expansion for x — > +00 [42] 

2 F 3 {a, b; c, d, e; -x) ~ Ax~ a + Bx~ b + Cx A/2 cos (2y/x + 

x (l + O^- 1 )) (14) 
with A = a + b — c — d — e+| and where the constants A, B, C are given by 
T(c)T(d)T(e) T(b-a) 



A = 
B = 

c = v^rr(a)r(6) 

Inserting this into the expression for S, we find 



r(6) F(c - a)F(d - a)F(e - a) 

r(c)r(rf)r(e) r(a - 6) 

T(a) r(c-6)r(d-6)r(e-6) 
T(c)T(d)T(e) 



(15) 



S ~ A (w*) a '"° + 5 (w*) ' -0-1 + (7 M^ 7 * sin (ut + | (a - A*/*)) (16) 

where the constants A, 5, (7 are proportional to A, B, C. The second term in (16) is 
always non-leading. The other two terms, however, will determine the functional form 
of the scaling function for y = ut sufficiently large. The treatment of x' is analogous. 

We can summarize the content of this section by listing the scaling functions which 
occur in (8), together with their leading behaviour as y — > oo 



X 'M = /„5(l-a',^-a) y 



l-a 



X 2^3 



1 - a' 2 - a' _ 3 1 - a - a' \ R 2 

2 ' 2 ' 2' 2 + 2z' _ 



a - a' Ar 

2 " 2z' 4 



(17) 



7T 



COS 



7ra 



r(i + a') 



-1 



+ / r 



-a y 



sin I y + 



7T 



[a - \ R /z]j 



XM = foB (-a', ^ - a) j,- 



-a' 1 — a' 1 —a — a' A# 1 — a — a' A# y 



7T 



sin 



2 



2z 



r(i + a') 



2 + 2z' 4 
'A 



(18) 



+ /o r ( ^ - a j y"W* cos (y + | [a - \ R / z\) 



We see that there appear terms which decrease monotonously with y but that 
there are also oscillating terms. They are described by different exponents and must be 
extracted by a different experimental setup. The oscillating terms follow the oscillations 
of the external field and the decrease of the oscillation amplitude gives a direct access 
to the exponent \r/z. On the other hand, in many experiments the data are averaged 
over one or several periods of the external field. For y sufficiently large, the contribution 
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of the oscillating term in eqs. (17,18) vanishes after averaging and then only a simple 
algebraic component remains, which permits to extract the exponent a — a'. For period- 
averaged data or else if Xr/z > a — a', the leading behaviour for large arguments is 

M ~ xSUO ~ y a '~ a (is) 

and the scaling relations (9) follow. 



3. Tests 

3.1. Numerical simulations 

We now compare the theory of the alternating susceptibility developed in the previous 
section with numerical simulations performed on the critical three-dimensional Ising 
spin glass. The Hamiltonian of the model is given by 

^ = -J2 J ^ a ^ ( 2 °) 

(i,3) 

where the nearest-neighbour couplings Jij are random variables taken from a binary 
distribution, i.e. Jij = ±1 with equal probability. The classical Ising spins Oi = ±1 
characterize the local magnetization at the sites i of a simple cubic lattice. In the past, 
numerical studies [43, 44] investigated the static critical properties of this model which 
undergoes a continuous phase transition at the temperature T c 1.19 (setting fc B = 1). 
Recently, we studied the ageing behaviour of this critical spin glass [13]. There we found 
clear evidence of a power-law scaling in the two-time correlation and integrated response 
functions. Here we shall need the exponent estimates a = 0.060(4) and X R / z = 0.38(2) 
obtained from a scaling analysis of the thermoremanent magnetization. 

In order to study the alternating susceptibility far from equilibrium we prepared the 
system in an uncorrelated initial state (corresponding to an infinite initial temperature) 
before quenching it to the critical point at time t — 0. At the same time an external 
oscillating magnetic field 

h(t) = h cos out (21) 

was switched on, with it's amplitude fixed at ho = 0.05 which is well inside the linear- 
response regime. We consider different values of the angular frequency uu = 2n/p with 
p ranging from 50 to 1600. Typically, systems containing 50 3 spins were simulated. 

Numerically, the in-phase and the out-of-phase susceptibilities are given by the 
expressions [45] 

x"(uj,t) = m(t) smut 
x'(ou,t) = m(t) cos cut 
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Figure 1. Scaling of the dissipative part ^'(wi) of the alternating susceptibility 
as function of the scaling variable Lot for different angular frequencies lj — 2n/p with 
p = 1600, 800, and 400. The full curve is the theoretical prediction (17) with f = 0.002 
and a' = —0.70 but which has also been shifted horizontally by y — > y + Ay, with 
Ay = —0.45, see text. Statistical error bars are smaller than the symbol sizes. 



with m(t) = 111 or der to access the scaling parts of x" an d of x' i we m ust 

i 

first subtract the equilibrium parts x'l an d Xx- We therefore carried out longer runs 
where we let the system equilibrate for typically a few ten thousand time steps before 
switching on the oscillating field. When no change in the amplitudes of the alternating 
susceptibilities were observed, we identified these data with the equilibrium parts x'l an d 
Xi- The ageing parts of x" an d of x' discussed in the following result from averaging 
over 2500 different runs with different couplings, different initial states and different 
realizations of the thermal noise. 

In figures 1 and 2 we test the expected scaling behaviour of the ageing part, see 



For the larger values of p, corresponding to the smaller values of u, we observe a very 
good data collapse for both quantities, which furnishes clear evidence in favour of a 
power-law scaling at T = T c . Here we used the value a = 0.060(4) determined previously 
from the decay of the thermoremanent magnetization [13]. For smaller values of p, the 
collapse is less good, which presumably means that for the corresponding values of u 
the dynamical scaling regime is not yet reached. 

The data shown in these two figures can in principle be compared directly with 
the analytical predictions (17) and (18). We point out, however, that the positions of 



eq. (7) 



X2 = and X2 = X^t). 



(22) 
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Figure 2. The same as in Figure 1, but now for the dispersive part x^O*^)- The full 
curve is the shifted theoretical prediction (18) with fo = 0.002 and a' = —0.70. 



maxima of x'2 an d X2 are shifted by the amount Ay « —0.45 when compared with the 
positions obtained in the analytical treatment of the previous section. The origin of 
this small shift is not completely clear to us. It is possible, however, that the cross-over 
between the equilibrium and the ageing regimes is not almost instantaneous in contrast 
to what we assumed in the derivation of eq. (10). In any case, we shifted the scaling 
variable in the analytical curves correspondingly, in order to make a comparison between 
analytical prediction and numerical data possible. 

As both a and Xr/z are known from our earlier investigation [13] the only free 
parameters in this comparison are the amplitude fo and the exponent a'. Our final 
estimates are 



and the fit is compared to the data in the figures. These values (23) of the parameters 
describe consistently both x'2 an d X2- Although some discrepancies are observed for 
small values of y = out, the overall agreement between the simulation and the shifted 
theoretical prediction is very good. The numerical data therefore support the scaling 
approach presented in the previous section. 

3.2. Comparison with experiments 

We now turn to a discussion of existing experimental results on the scaling of the 
alternating susceptibility. 

A detailed discussion of a possible scaling of the alternating susceptibility was 



a' = -0.70(3) 




(23) 
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presented by Suzuki and Suzuki for the short-ranged Ising spin glass Cuo.5Coo.5Cl2- 
FeCl3 graphite bi-intercalation compound [37]. After a rapid quench to below the glass 
temperature T g = 3.92(11)K, they measure x"(w,t) for fixed u> and find that their 
(period-averaged) data are well-fitted by the power-law 

FM) = xSM + ^»r'' (24) 

While the fitted exponent b" depends only slightly on u, they further show evidence for 
a power-law A"{u) = AqU~^" and find that ". . . the value of /i" is almost the same as 
that ofb"" [37]. In this way, they arrive at the cut-scaling form 

r(uJ,t) = X » + A>UuJt)- b " (25) 

and a similar form for x'(w,t) where ". . .b' and b" are of the same order at the same 
temperature" [37]. Experimentally measured values of the exponent b" (and also b') of 
some materials are collected in table 1. We now compare the experimental results of 
Suzuki and Suzuki [37] with the theoretical scaling form (8). First, the experimental 
evidence for a pure u;t-scaling indicates that the exponent a must indeed be very small. 
Second, when considering the leading behaviour for y = ojt large (their data go up to 
y < 10 6 [37]) and recalling that the experimental data are averaged over at least one 
period of the external field, we can read off 

b" = a-a' , n" = -a' (26) 

and the observed [37] near equality b" « fj," is again consistent with a being negligibly 
small. Third, the available data are consistent with the theoretically requested relation 
b' = b" . Forth, the reason for the experimentally observed sudden jump in b' and b" for 
smaller T is not yet understood but we mention that a similar phenomenon also occurs 
in certain relaxor ferroelectrics [40, 46]. 

Similar values of b" were observed for several other materials, quite independently of 
the precise physical nature as can be seen from table 1, but the errors are still too large 
to permit a discussion of the universality of the exponents. However, the experimental 
data are in many of these systems at least as well described by a logarithmic scaling 
as expected from the droplet theory [48, 32, 49]. Furthermore, in several systems also 
strong deviations from a simple u;t-scaling were found, see [46]. Finally, we mention 
that in systems like /3-hydroquinol-clathrate [41] or even simple liquids like glycerol [50] 
a power-law dependence of the form x'l gc ~ t a or xl gc ~ ^ a was observed. All in all, it 
is not yet completely understood what precise conditions are needed such that a clear 
power-law scaling in ut can be observed. 

Lastly, we see that the value b" ~ 0.7 obtained from the critical 3D Ising spin glass 
with binary disorder is very far from the values b" w0.1- 0.2 found experimentally. 
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Table 1. Measured values of the exponents b" and b' in several glassy materials, using 
the scaling form (7). Here T g stands for the glass transition temperature and T is the 
temperature where the data were taken. For Feo.sMno.sTiOa and CdCr1.7Ino.3S4 the 
relation b' = b" was assumed. 



Material 


T g [K] T[K] b" b' 


Ref. 


Cu . 5 Coo.5Cl 2 -FeCl 3 
- GBIC 


3.92(11) 3.25 0.01(4) 0.08(3) 
3.5 0.017(32) 0.05(2) 
3.75 0.16(3) 0.20(2) 
3.85 0.15(3) 
3.95 0.16(4) 0.20(2) 


[37] Ising spin glass 


Fe . 5 Mn . 5 TiO 3 


20.7 15 0.14(3) 
19 0.14(3) 


[49, 48] Ising spin glass 


CdCrx.7Irio.3S4 


16.7 12 0.18(3) 
14 0.18(3) 


[49, 48] Heisenberg spin glass 


CdCr 2 Jn 2 _ 2a; S4 x = 0.95 
x = 0.90 


70 8 0.2 
67 0.2 
50 42 0.20 


[51, 48] disordered ferromagnet 
[47, 48] 


Pb(Mg 1/3 Nb 2/3 )0 3 


~ 220 < 220 0.17 


[40] relaxor ferroelectric 



4. Conclusions and discussion 

The objective of this work has been to investigate to what extent the ageing behaviour 
of the dispersive and dissipative parts of the alternating susceptibility in glassy systems 
may be described in terms of some simple ideas borrowed from the ageing of simple 
magnets without disorder. These are 

(i) the clear separation of the stationary and the ageing regimes, which goes into the 
derivation of eq. (10). 

(ii) the hypothesis of a single essential length scale L(t) ~ t x l z growing algebraically in 
time. 

(iii) the extension of this dynamical scaling to a local scale- invar iance which leads to 
the simple form eq. (4) for the two-time response function R(t, s). 

Taking these assumptions as working hypothesis, our results are as follows: 

(i) The often-studied simple wt-scaling (7) in x" an d x' should be slightly generalized 
to the scaling forms (8) or equivalently 

x "(u;,t) = x» + t- a x2(ut) with xi(y) = y a X2(y) (27) 

and similarly for x' ■ The smallness of the exponent a makes it conceivable that 
the slight 'subageing' found in many experiments might be taken into account this 
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way.* Tests of this idea in real materials would be welcome. 

(ii) The exponents b' and b" defined from X2(y) ~ V~ h ' an d x'iiv) ~ V~ h " satisfy the 
scaling relation 

b' = b" = a - a' (28) 

which hence relates exponents found in an oscillating field to those describing the 
two-time response R(t,s). Since a' appears to be a new, independent exponent, 
this suggests that the exponent b' = b" should be independent of the exponents 
a,b, Xc/ z, Xr/z habitually used in describing ageing. This question could be 
addressed, independently of the hypothesis of local scale-invariance and the form 
(4), through field-theory methods as developped in [52]. 

(iii) When testing the explicit dispersive and dissipative scaling functions (17) and (18) 
on the 3D critical Ising spin glass with binary disorder, we found that the form of 
both is reproduced by the data for the following exponents 

a = 0.060(4) , a' = -0.70(3) , ^ = 0.38(2) (29) 

(and with a normalization constant f ~ 0.002) where we took over the values of 
a and X R / z from our earlier analysis of R(t,s) in this model [13]. We stress that 
this agreement is only found if the curves (17) and (18) are shifted y — > y + Ay 
by a constant Ay ~ —0.45. Even then, for small values of y the simulated scaling 
functions show a small but systematic deviation from the theoretical prediction 
which remains to be understood. 

The origin of this shift is unknown to us, but it might be due to a rather slow 
cross-over between the stationary and the ageing regimes which is not captured by 
eq. (10). Understanding this point is an important open problem. 
We think it is remarkable that the ageing dynamics of a system as complicated 
as a critical spin glass can be captured by our relatively simple hypothesis of a 
dynamical symmetry. 

(iv) The value b" ~ 0.7 we found in the 3D Ising spin glass (with binary disorder) is far 
from the values obtained in experiments, see table l.Jj If the present result should be 
confirmed, it might furnish a further hint towards a fundamental difference between 
simple spin glass models and experimentally realized glassy systems [36]. 

On the other hand, the theoretically predicted relation b' = b" is fully consistent 
with the available experimental results. 

Acknowledgements: We thank G.F. Mazenko for helpful correspondence and 
O. Petracic for a useful discussion. This work was supported by the Bayerisch- 
Franzosisches Hochschulzentrum (BFHZ) and by CINES Montpellier (projet pmn2095). 

* Working with (27) rather than with (8) has the advantage that for y 3> 1 the first term in the 
asymptotic expansion only depends on the exponent a', while the second one only depends on a — Xr/z. 
jj Our previous investigation of the thermoremanent magnetization [13] did not permit us to determine 
the value of a'. 
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Appendix A. The autoresponse function and extended local scale-invariance 

Local scale-transformations are extensions of the dynamical scale-transformation r i— > 
br, t i— > b z t of space-time towards variable rescaling factors b = b(t, r) such that 
conformal transformation in time 1 1— > (at + f3) / {pft + 5) with a5 — j3 r y = l are maintained 
[14]. For any given value of z, such infinitesimal transformations have been explicitly 
constructed and shown to furnish a dynamical symmetry of the generalized diffusion 
equation (d t + d*) ip(t, r) = 0. Here we are interested in applications to ageing behaviour 
of the autoresponse function 

R(t,s) iM 



5h(s, r) 



0(t,r)0(s,r)) (Al) 

h=Q 



where is the order parameter, h the conjugate magnetic field and the associate 
response field in the context of the Martin-Siggia-Rose theory, see e.g. [53]. Therefore, 
we have to restrict to the subalgebra with time-translations excluded [14, 15]. Since 
space translations are included in the set of local scale-transformations, the scaling of 
R(t, s) will be entirely given in terms of merely two generators which we write as 

X = -tdt - - , X 1 = -t 2 d t --(x + Ot (A2) 

z z 

where x is the scaling dimensions of the field on which the generators X n act. We have 
observed earlier [17] in the special case z = 2 that a further constant £ can be introduced 
without changing the commutator relations of the subalgebra of local-scale invariance 
under consideration and now extend this to arbitrary values of z. Previous treatments 
of the question had admitted £ = £ = from the outset. The covariance of R(t, s) is 
now expressed as usual [14, 15] 

X R = (-td t ~ sds ~^~-^) R & s ) = 



XiR = (-t 2 d t - s 2 d s -^(x + i)t-- z (x + t) s^j R(t, s) = 



(A3) 



where x and £ refer to the response field 0. To solve these, change variables into u = t — s 
and v = t/s. Then, with R(t, s) = R(u, v) 

x ~\~ x \ — 
ud u H — j R(u, v) 

, v r - .r ■ 2i 1 x - x + 2£ . - 

« 4+ 7 - + 7 -\R(u,v) = A4 

v — 1 z v — 1 z 




Scaling and ageing of the alternating susceptibility in spin glasses 15 

The solution to these is found in factorized form R(u, v) = f(u)g(v) and we find, after 
having returned to the variables t and s 

/. \-l-a' / ,\l+a'-\ R /z 

R(t, s) = rvs- 1 -* - - 1 - (A5) 



where 



x + x 
1 + a = 

z 



i + a , = x + S+ 2Z + 2Z 

z 

\ R = 2 (x + £) 

which is the form (4) stated in the text. We finally observe that if £ + £ = 0, we recover 
indeed a = a'. 



Appendix B. On the Ohta-Jasnow-Kawasaki approximation 

We briefly present the derivation of the scaling of the response function R(t, s) in a simple 
analytically tractable scheme which is close in spirit to the Ohta-Jasnow-Kawasaki 
(OJK) approximation. Following the ideas of Berthier, Barrat and Kurchan [21] and of 
Mazenko [22] which in turn are based on a calculation by Bray [54], one considers the 
integrated response in the zero-field cooled protocol. When perturbing the system by 
a spatially random field h of magnitude ho, the zero-field-cooled (ZFC) susceptibility 
reads [21, 54] 

(h(r)<f>(t,r)) [2(h(r)m(t,r)) 



X(t, s) = V V T 2 V ' " — \ ~ ^ /rr^r (Bl) 



m 2\ 



where (j) is the order parameter and it is assumed that for late times, one can approximate 
4> ~ sign (m). At zero temperature, the auxiliary field m should satisfy the equation of 
motion 

— - = V 2 m - n a n b V a V b m + h\Vm\. (B2) 
at 

In the spirit of the OJK-approximation, one makes [54, 2] the simplifications n a n b — > 
5ab/d (circular average) and \Vm\ — > ((Vm) 2 ) 1 / 2 . The equation of motion then becomes 
[21, 54] d t m = DV 2 m + h((Vm) 2 } 1/2 with D = (d - l)/d and assuming the fields m 
and h to be gaussian, it is found that [21, eqns. (19,20)] 



which is the starting point of our analysis. 
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Performing the integration in fc-space, one obtains [22] 

X (t, s) = cste. J\u u^ d+1 ^ 2 (^j ' (i- - ij ' = J\u R(t, u). (B4) 

This integral becomes singular near the upper limit u ~ t and we shall reconsider this 
below. Before we shall do this, we read off the autoresponse function 

R(t, u) = u"^/ 2 f om {t/u) , f OJK (y) = / y dt \y - iy d/2 (B5) 

from which we easily recover, using eq. (4), the values of the exponents a, a' and \r in 
the OJK-approximation (where z = 2) as stated in eq. (5). Eq. (B5) is also recovered 
in the gaussian theory of phase-ordering [23]. 

A different conclusion was reached by Berthier et al. [21] albeit starting from the 
same eq. (B3). They quote x(t, s) ~ s~ 1 / 2 F(t/s) for d > 2 which would mean a = 1/2, 
provided of course that the naive scaling law x(t, s) = s~ a f x (t/s) could be used. We 
must therefore reconsider the singularity in eq. (B4) for x(^ s )- Following [21, 22], one 
introduces a cut-off parameter A 2 (which should be sent to zero at the end) and writes 
instead of (B3) 

Performing first the integral over k and changing variables, this becomes, up to a 
normalization constant 

X (t, s) ~ t^' 2 f dv v-^ 2 ^ (l-v + A 2 /t) ~ m (B7) 

J s/t 

and we must analyze the contribution of the integrand near v ~ 1. We decompose the 
domain of integration J g , f = J s ^ £ + J l £ . In the first term, we let A 2 — > 0, and v ~ 1 in 
the first factor of the second term. Then, for d > 2 

X (t, s) ~ t-^' 2 f^dv v-^/\l - v)- d / 2 

J s/t 

+ r( d - 1 )/ 2 f dv (i-v + A 2 /ty d/2 

Jl-e 

= t- a F(s/t) + dr 1 / 2 (A 2 ) 1 "" 72 + c 2 r a (e + A 2 /t) l - d/2 

= s' 1 / 2 ■ xoo (t/sy 1/2 + s- a f x (t/s) (B8) 

where the value a = (d — l)/2 was used. Here Ci,2 and Xoo are constants and F and f x 
are scaling functions. In this way, not only we recover the leading term already found in 
[21] but we also see that the contribution coming form formally integrating the scaling 
behaviour of R(t, s) merely gives rise to a sub-leading correction. A similar argument 
can be applied to the case d — 2 and produces logarithmic corrections. 
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Consequently, the mere observation of a scaling law for x(t, s) in the scaling regime 
t, s — > oo with y = t/s fixed is not enough to reliably extract the exponent a. Indeed, the 
term expected from naive scaling s) ~ s~ a merely arises as a short-time correction 
to the leading long-time behaviour of s). It has already been pointed out in the 
context of the Ising model quenched to T < T c that a similar dominant term not simply 
related to the autoresponse exponent a occurs in the ZFC-susceptibility for d < 3 [55]. 
The straightforward use of x(t, s) may hence lead to erroneous values of the exponent 
a and it is safer to avoid using x(t, s) altogether. 
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